Exact solution of the spin-isospin proton-neutron pairing Hamiltonian 
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The exact solution of proton-neutron isoscalar-isovector (T=0,1) pairing Hamiltonian with non- 
degenerate single-particle orbits and equal pairing strengths (<7t=i=St=o) is presented for the first 
time. The Hamiltonian is a particular case of a family of integrable SO (8) Richardson-Gaudin 
(RG) models. The exact solution of the T=0,1 pairing Hamiltonian is reduced to a problem of 4 
sets of coupled non linear equations that determine the spectral parameters of the complete set of 
eigenstates. The microscopic structure of individual eigenstates is analyzed in terms of evolution 
of the spectral parameters in the complex plane for system of A=80 nucleons. The spectroscopic 
trends of the exact solutions are discussed in terms of generalized rotations in isospace. 
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The exactly solvable models introduced by Richardson 
[H and by Gaudin [3] belong nowadays to classic theo- 
retical tools in mesoscopic physics. Indeed, these models 
based on the rank 1 SU(2) algebra for fermions or the 
SU(1,1) algebra for bosons were applied to a large variety 
of quantum many-body systems including the atomic nu- 
cleus, superconducting grains, cold atomic gases, etc., see 
review article Q and refs. therein. Recently, we have ex- 
tended the Richardson-Gaudin (RG) models to the rank 
2 algebras: SO(5) (isovector pairing Q), SO(3,2) (F spin 
1 boson pairing 0]), and SU(3) (interacting three level 
atoms @). 

In this letter we will derive for the first time the exact 
solution for the rank 4 SO (8) RG integrable model with 
non-degenerate single-particle (sp) spectrum and arbi- 
trary degeneracies. As a particular realization of the rank 
4 SO (8) RG model we will consider the nuclear isoscalar- 
isovector (T=0,1) pairing Hamiltonian introduced for a 
single degenerate shell in Ref. Q and further developed 
in [1, H| . We will solve the model for a realistic case of 
A— 80 nucleons moving in fourfold degenerated equidis- 
tant sp spectrum working in T—0,l pair representation 
of the SO (8) algebra. It should be mentioned that other 
representations like the Ginnocchio model [l(| can lead 
to interesting exactly solvable models in nuclear structure 
as well as to models of spin 3/2 cold atoms [Tl| . 

The study of proton-neutron (p-n) pairing has gained 
a renewed interest due to the new generation of 
radioactive-beam facilities that will open the access to 
proton-rich nuclei close to the N=Z line. In spite of 
vigorous activity in this field, see [l2| and refs. therein, 
the fundamental questions concerning the basic building 
blocks and experimental fingerprints of the p-n pairing 
are still a matter of debate. So are the theoretical prob- 
lems concerning generalization of well established nuclear 
pairing models to include p-n pairing, proper treatment 
of isospin degree of freedom or a-like clustering. All these 
problems set clear motivation for realistic exact-model 
studies of the p-n pairing undertaken in this work. 

Let us begin our derivation by introducing the 28 gen- 
erators of the SO(8) algebra [tJ: three (T=l,5=0) and 
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P Tl = (Ptjt, and D CTi = 

(D^)^, where the triads in the couplings represent, re- 
spectively, angular momentum, isospin and spin. The 
fermionic operators a|. mro . create a particle in the orbit 
li with projection to, isospin r and spin a. The SO(8) 
algebra is completed by the 16 particle-hole operators: 
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These 16 operators 



close an U(4) subalgebra of SO (8) and include the num- 
ber operators for the four different nucleon types in the 

Orbit li. ]M Tcr j — ^Ttx,T<r;i • 

The general procedure for solving RG models for ar- 
bitrary simple Lie algebras has been developed in refer- 
ences [III, [3| ■ For each algebra it is possible to derive a 
set of L quadratic integrals of motion defining the inte- 
grable model. It is also possible to derive the complete 
set of common eigenstates and eigenvalues, which con- 
stitute the exact solution of the model. Here L is the 
number of copies of the algebra that we associate with 
the number of orbits U. For simplicity we will be con- 
cerned here with a particular Hamiltonian, the T = 0, 1 
proton-neutron pairing Hamiltonian that arises as a par- 
ticular linear combinations of the integrals of motion of 
the SO (8) RG model: 
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where Nj is the number operator of the orbit li. 

We would like to emphasize here that the Hamilto- 
nian (JTJ) has equal strength for T=l and T—0 pairing. 
As a consequence, there is a conserved U(4) symmetry 
defined by the generators C Tl<Tlj7 - 2Cr2 = £\ C Tl<TljT2Cr2:l . 
Therefore, the eigenstates are organized in degenerated 
U(4) Wigner multiplets. For a given number of nu- 
cleons A, these U(4) multiplets can be classified using 
Young tableaux. Each multiplet is defined by a par- 
tition of A in 4 numbers, [ft.1/12^3^4], constrained by: 
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J2ii 2 h + 1) > hi > h 2 > h 3 > hi > 0. The labels hi 
are related to the number of particles in the total U(4) 
Lowest Weight State (LWS). For instance, if we relabel 
the U(4) operators according to the rule 1 = n|, 2 = n\, 
3 = pi and 4 = pi, the U(4) LWS can be defined as the 
state which satisfies: C a ^\LWS) = 0, Va < [3. For this 
choice of LWS, the corresponding U(4) Young tableau 
is given by: [N n iN n ^N P iN p ^\. The spin and isospin 
of this LWS are simply 2S = N nl +N nl ~N nl -N nl and 
2T=N nL +N nr -N pl -N pV 

As stated above, the eigenvalues of Hamiltonian (fl]) 
can be derived from the exact solution of the SO (8) RG 
model. They are: 
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where Ui is the seniority of level i, i.e., the number of par- 
ticles in level i not coupled in T—l or T—0 pairs. The 
parameters e a satisfy the generalized Richardson equa- 
tions: 
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where [/ii;i^2;i^3;i^4;i] is the Young tableau of the re- 
duced U(4) irrep defined by the unpaired particles in 
the i-th orbit. These U(4) labels are constrained by 
the conditions: 2^ + 1 > hi-i > h 2 -i > h^-i > hi-,% > 
0. In terms of these labels the seniority of level i is 
Ui = J2khk-,i- The rank of the RG models defines the 
number of different sets of spectral parameters. SO (8) 
is a rank 4 algebra, hence there are 4 sets of spec- 
tral parameters. The number of spectral parameters in 
each set is determined by the reduced labels and those 
of the total U(4) Wigner multiplet: Mi = (A - u)/2, 
M 2 = h 3 + h 4: - J2i( h 3;i + ^4;i)> M 3 =hi-Y,i hi;U and 
M 4 = (A- 2h 1 )/2 - (u - Ei 2fci ;i )/2, with u = '£ i Ui . 

The first set of spectral parameters comprises the usual 
pair energies e Q of the SO(8) algebra. The other three 
sets, composed by the spectral parameters cu a , rj a and 
7 Q , are associated with the U(4) subalgebra of SO (8). 



While the eigenvalues depend only on the parameters e a , 
the corresponding eigenfunctions are determined by the 
parameters of the four sets. The complete set of solutions 
of the Richardson equations defines a basis which spans 
completely the Hilbert space of sates with the same U(4) 
Wigner quantum numbers [/iift.2/2.3/14]. 

Even though the set of non linear coupled equations 
([3]) seems to be extremely complex, we will show how 
it is possible to obtain numerical solutions within a non 
trivial example of A— An nucleons with n a positive inte- 
ger moving in a set of non-degenerate I = single particle 
orbits. Other cases could be handle following a similar 
procedure. Before describing the numerical strategy, it 
will be useful to consider the lowest energy LWS config- 
urations for even and odd isospin in the g=0 limit. 

For even T the lowest (A ~ 2T) /4 levels are filled with 
4 nucleons, and the following T levels with a pair of neu- 
trons. All particles are paired and the corresponding se- 
niority quantum numbers are 0. In the case of T odd, the 
levels (A-2T+2)/4and (A+2T+2)/4 have one unpaired 
nucleon (u( j4 -2T+2)/4 = U(A+2T+2)/4 = !)• This state 
can be considered a p-h excitation that evolves to a two 
quasi-particle state in the superconducting phase. As a 
consequence, the number of pair energies e Q is Mi = A/2 
in the even T case and Mi = A/2 — 1 in the odd T case. 
In this limit the pair energies take the values e Q = 2e; 
according to the same pattern. 

In the weak coupling limit (j « 1) the Richardson 
equations ([3]) decouple into independent sets of equa- 
tions, each one related to the single particle level partially 
or fully occupied in the g — limit. These equations can 
be solved analytically. The 4 sets of spectral parame- 
ters obtained in this way are used as initial guess for an 
iterative procedure in which the coupling constant g is 
increased step by step using the previous solution as the 
initial the guess. 

As is well known, the main obstacle in solving the 
Richardson equations even in the SU(2) case, is the ap- 
pearance of singularities at some critical values of the 
pairing strength due to crossings in the real axis of single 
particle energies and spectral parameters [H, . This 
problem is even worst in the SO (8) model having 4 sets 
of spectral parameters. In order to avoid these numerical 
instabilities we introduce an alternate imaginary term in 
the sp-energies ej — > + (— ) l Ay / — T, which breaks the 
time reversal symmetry and moves the solutions of ([3j) 
away from the real axis [171 ] . The system is then evolved 
from the initial guess at g <C 1 to the desired value of 
g. At this point we begin a second iterative process to 
set the imaginary term in the single particle energies to 
zero (A — ■> 0). This recursive procedure proved to be 
very efficient in solving the Richardson equations for ar- 
bitrary values of the coupling constant and for all the 
states considered. 

We will now demonstrate the ability of our procedure 
for solving large scale complex physical problems by pre- 
senting a numerical example for a system of yl=80 nu- 
cleons described by the Hamiltonian (jTJ) and moving in 
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FIG. 1: Complex plane representation of the pair energies e 
and wave function structure parameters w, 77 and 7, for the 
lowest-energy states with isospin T=0,4,8. The values of g 
are indicated on the top of the columns. Horizontal lines 
represent the sp-energies 2e; = (i — 1). Clusters representing 
alpha quartets (T=0) and collective p-p pairs (T=4 and T=8) 
are indicated by full and dashed ellipses respectively. 



a set of T=50 equidistant and fourfold degenerate levels 
(e* = (i- l)/2, Z 4 = with i = 1, ...,50). The advan- 
tage of the method is that it provides not only exact 
eigenstates (or spectroscopic information) but also allows 
for intuitive, pictorial representation of the microscopic 
structure of individual eigenstates. In the following we 
will briefly present both aspects of our model. 

Let us start with a brief discussion of microscopic as- 
pects of our solutions. Since the exact eigenstates are 
fully determined by the spectral parameters (e, u>,r], 7) 
the evolution of these parameters in the complex plane 
allows for tracking structural changes of individual eigen- 
states as a function of model parameters. This is demon- 
strated in Fig. [JJ where we show the spectral parameters 
in the complex plane for three different isospin states and 
two values of the coupling constant representing weak 
g = 0.16 and intermediate g — 0.22 pairing. As in the 
SU(2) RG model, the expansion of the pair energies e 
into the complex plane indicates the formation of corre- 
lated Cooper pairs [l8|. Simple counting shows that, for 
the two g values considered here, the overall number of 
correlated pairs is, irrespective of isospin, about 6 (15%) 
and 14 (35%) respectively. The rest of the pairs, still 
attached to the sp energies, remain almost uncorrelated. 

The expansion of the spectral parameters u>, rj and 7 in 



the complex plane follows the behavior of the pair ener- 
gies e, namely, they form parallel arcs to those of the pair 
energies and they arrange themselves into various cluster- 
like structures, see Fig. [TJ According to our choice of the 
LWS an isolated complex parameter e represents a col- 
lective T=l n-n Cooper pair; a cluster of one e, two w's, 
one 77, and one 7 represents a collective T=l Cooper p-p 
pair, while a cluster of two e's, two u'a, one 77, and one 
7 represents a correlated alpha-like quartet. This allows 
for an unambiguous interpretation of Fig. [TJ The exis- 
tence of correlated quartets is clearly visible in the T = 
panels of the figure for strong enough g. With increasing 
T, these clusters break apart and the net separation of 
the arcs of the pair energies increases. Moreover one of 
the arcs is formed by isolated pair energies e while the 
other one is constituted by pair energies forming clus- 
ters with two w's, one 77, and one 7 spectral parameters. 
Physically, it implies a quenching of isoscalar pairing and 
the formation of two separate conventional n-n and p-p 
pairing condensates. 

Let us turn now to the spectroscopic consequences 
of the observed microscopic processes along the nuclear 
symmetry energy (NSE) curve E(T). In the following we 
will show that these processes form a systematic pattern 
which can be nicely interpreted in terms of generalized 
rotations in isospace in the spirit of iso-cranking model 
of Refs. [lj| . Let us recall that according to that model 
the NSE splits into two structurally different even-T and 
odd-T branches of an iso-rotational band which can be 
conveniently parametrized as: E (e \T) = T(T + X)/2J T 
and £(°)(T) = T(T + \)/2J T + AE exc respectively Here 
Jt stands for moment of inertia in the isospace (iso-Mol) 
while A determines strength of the linear term ~ T which 
is often called the Wigner energy. The odd-T sequence 
is shifted up with respect to the even-T branch by a 
two-quasiparticle [2qp] excitation energy AE exc . Note 
a beautiful analogy to the spatial collective rotation in 
even-even nuclei where odd-spin branch is also built upon 
2qp excitation. 

The calculated inverse of the iso-Mol (1/Jr) versus 
T, the primary characteristic of the iso-rotational mo- 
tion, is shown in Fig. [2^. Apart from the SO (8) solution 
also two limiting SO (5) cases invoking only isoscalar and 
only isovector pairing are depicted [4[. Note that in ac- 
cordance to the iso-cranking model: (i) all curves con- 
verge to the sp splitting 1/J T -> Se = 1/2 @, HH (ii) 
the T=l paring represents almost perfectly rigid rotation 
with 1/J T w Se irrespectively on T [HHJ (Hi) the SO(8) 
and T=0 pairing curves show a characteristic reduction 
of the Jt at low-T due to isoscalar pairing collectivity 
similar to the well recognized reduction of the spatial 
Mol caused by isovector superfluidity. The increase of 
Jt versus T reflects disappearance of the isoscalar pair- 
ing collectivity caused by fast iso-rotation which tends to 
recouple isoscalar (anti-parallel coupled isospins) pairs in 
analogy to the well known Coriolis anti-pairing effect. 

The quantity depicted in Fig. [2]b is directly related to 
average difference in pairing correlation energy E corr = 
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FIG. 2: (a) Inverse of the iso-Mol (1/J-r), (b) signature split- 
ting Sex and (c) linear term enhancement factor A versus T 
for pure T=0 model (circles), pure T=l model (diamonds) 
and the SO(8) model (triangles). Filled (open) symbols re- 
fer to even-(odd-)T branches of E(T). The calculations were 
done for g = 0.16 except for gray triangles in the lowest panel 
which mark the SO(8) solution for g — 0.22. 



(H) — E(T) (() denotes the expectation value in the g=0 
limit) between even and odd-T branches. This quantity, 
known also as signature splitting, is defined at odd-T as: 
Se T ee (E^(T+l) + E^(T~l))/2-E ( ^(T). At high- 
T, where 1/Jt ~ Se, signature splitting equals Sex — 
~AE exc + Se/2 irrespective of A. In the case of pure 
T=0 pairing building up the isospin proceeds through 
the isoscalar pair breaking. The correlation energy drops 
down and the solution goes smoothly over to the sp limit 
where AE exc — > Se/2. Consequently Sex — > as shown 
in Fig. [2b- In the T=l pairing case the signature splitting 
is almost constant and equal Sex ~ —1.5. In this case we 
deal with rigid iso-rotation and odd and even-T branches 
are shifted by a constant energy reflecting a difference 



between correlation energies in odd (seniority two) and 
even-T (seniority zero) states (2qp energy). Finally, the 
SO (8) curve goes smoothly over to the T=l case as the 
isoscalar pairing disappears with increasing T. 

Fig. [2]; shows the linear enhancement factor calculated 
as A ee 2E T e ' Jt/T — T where Jx stands for mean value of 
the iso-Mol. While the T=l pairing yields A w 1, strong 
enhancement of the Wigner term due to the isoscalar 
pairing is clearly seen as anticipated [22]. In the SO(8) 
case A reaches the Wigner supermultiplet limit A ~ 4 for 
large g (gray triangles) and drops with decreasing g as 
well as with increasing T reaching unity for large T. 

In summary, we have presented the exact solution of 
the RG model associated to the SO (8) algebra in the 
context of nuclear n-p pairing with equal strength for 
the T=l and T=0 interaction components. We have 
briefly discussed a new technique for solving the Richard- 
son equations which has a potential to become an in- 
valuable tool in studying integrable models for binary 
mesoscopic systems. The first application to the nuclear 
n-p pairing is discussed from both microscopic as well as 
spectroscopic points of view. In particular, it is shown 
that T = wave function shows alpha-like quartet struc- 
tures that can be recognized by the formation of clusters 
of spectral parameters containing two pair energies. At 
high T these alpha-clusters dissolve and two separate p- 
p and n-n superfluid condensates are formed. Spectro- 
scopic consequences of these microscopic processes are 
discussed and interpreted in terms of of generalized rota- 
tions in isospace. It is shown that the exact solutions 
follow nicely the general trends predicted by the iso- 
cranking model. 
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